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FISSION INDUCED PLASMAS 


By 

Wynford L. Harries^ 


1. INTRODUCTION 

The grant commenced on December 15, 1976 and runs for one year. 
This report covers the period through June 30, 1977. 

The purpose of the work was to investigate the possibility 
of creating a plasma from fission fragments, and if possible to 
utilize the energy of the particles to create population inver- 
sion that would lead to laser action. Eventually it is hoped 
that the same medium could be used lor both fissioning and lasing 
thus avoiding ineff iciences in converting one form of energy to 
the other. 

It was originally intended that a research associate be hired, 
and two applicants with backgrounds in nuclear pumped lasers were 
interviewed; however, neither accepted the position. Accordingly, 
after June 15 the writer undertook the investigation; but prior 
to that (because of the academic year and other grants), only 
seven percent time was charged on this grant. 

One of the central problems in understanding a fission induced 
plasma is to obtain a model of the electron behavior, and some 
calculations are presented in section 2. The purpose of this sec- 
tion is to provide a handy compendium of processes for reference. 
Some sections, especially 2.1.1. and 2.1.2. are largely due to 
John Wilson, but are included for convenience and because the 
numbers are needed later. The calculations are simple, and 
intended to get orders of magnitude only in many instances. 


Professor, Department of Physics and Geophysical Sciences, 
Old Dominion University, Norfolk, Virginia 23508. 


1 


Section 3 describes possible future work, especially experi- 
ments for checking the parameters in section 2. 

2. ELECTRONS IN A FISSION INDUCED PLASMA 

We consider the electrons to be in a background gas of He^ 
for this discussion, but the results will also apply for small 
admixtures of other gases. 

2.1. Production of Electrons 

The electrons are produced by the products of the action of 
thermal neutrons on He^ : 

n + He^ p(0.57 MeV) + T(0.19 MeV) 

Both proton and triton ionize. 

2.1.1. Number of Electrons Produced Per Second 

On the average, 34 eV are needed to produce one electron-ion 
pair, and the electron receives an average of 8 eV kinetic energy; 
thus, 20 percent goes in kinetic energy (ref. 1). Hence, the 
number of electrons produced per neutron, Q, is 

Q = [(0.57 + 0.19)MeV/(34 + 8)eV] =: 2 x 10“* electrons. 

If the fast burst reactor produces a flux, F, of 10^^ 
neutrons per cm^ per sec, and they are absorbed with an average 
range r = 6 cm, then the number of electrons produced per sec 
per cc is 

I « QNr”^ = 3 X 10^° el per cc per sec. 
prod 

A more accurate estimate is based on a calculation by Wilson 
(ref. 1) for the energy deposited per neutron per cm^ , in unit 
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volume, which is 0.3 MeV per cc per neutron per cm^ for p = 3 
atmospheres. Multiplying this by the flux gives the energy input 
per cc per sec, namely, 3 x 10^^ eV per cc per sec. Assume that 
at one atmosphere the energy deposited is one-third as much, or 
10^^ eV per cc per sec. The production rate of electrions is 
this number divided by (34 + 8)eV or 2 x 10^^ el per cc per sec. 
Both these numbers roughly agree, and so we assume 

j « 2 X 10^® el per cc per sec. 

/ prod 

2.1.2 Initial Electron Energy Distribution * 

The fission fragments ionize the He^ and release electrons. 
For high energy fragments whose velocity exceeds that of the 
orbital velocity of the electron, the binding energy of the 
electron can be neglected, and the collision regarded as 
Rutherford scattering. The differential cross-section subtended 
is proportional to E"j, where is the energy of the fission 

fragment, and the electrons are expected to be in an initial 
energy distribution function, fi(e) « where e is the 

energy of the electron. The maximum amount of energy handed 
over to the electron is given by e « (2m/M)E-, where m is 
the electron mass, and M the mass of the fission fragment. 

The initial distribution function for electrons created by a 
fragment of energy E^ ^ is shown in figure 1 (solid curve). 

The picture fails around 1.5 e^, where is the ionization 

energy, or around 30 eV, as there are no more excitation levels, 
and the curve goes to a maximum and falls off at lower energies. 

As the fission fragment slows down to an energy E. the 

X f * 

cutoff drops to E „, and the distribution is now as shown by 

O j A 

the dotted curve. The distribution function for all the electrons 
produced by one fission fragment is the sum of such curves, and 
will have a slope steeper than e”^ 


* This description is due to J. Wilson. 
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Figure 1. 


Initial distribution function fi(e) 
for electrons; the curves are not 
normalized . 



This is the energy distribution function of electrons just 
created by ionization; thereafter they lose their energy by 
inelastic and elastic collisions. 

2.2. Electron Energy Losses 

2.2.1. Energy Losses by Inelastic Collisions 

In the range 20 < e < 2000 eV, the most rapid energy loss rate 
for electrons is due to inelastic collisions of ionization and 
excitation of both neutrals and ions. The rate of change of the 
energy e of one electron with time is 

^ = _ViEi - V 2 E 2 - ... = -vE (1) 

where is the collision frequency for level n, and E^ is 

the energy given up by the electron. The quantity ^ is an 
average for all levels, where 

V 

V = -jp = N aVr (2) 


and v^ is the electron random velocity, X is the mean free 
path, N is the density of targets, and a the cross section 
for inelastic collisions. Here we are concerned with orders of 
magnitude only and assume X = constant. Then 


1/2 

v^ = (2£/ra) 


(3) 


and 


d£ _ 
dt 


No 




2 V E 

m/ X 


1/2 

e 


(4) 
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Thus the energy of the electrons varies with time according to 



where is the initial energy of the electron and is a 

characteristic time for inelastic energy loss, 


T. 

in 


o 

E 


( 6 ) 


The excitation cross section for He above 19 eV is 0.5 x 10“’ crn^ 

(ref. 2), and the number of atoms per sec at 760 torr is 2.7 x lo^®. 

Hence X = (Ncr)“^ = 7.4 x lO”^ m. The value of T^^^ depends on 

the choice of e , and if e w 1 KeV, and E ~ 20 eV, then 

o’ o 

T. « 4 X 10" sec. 
in 

In the energy regime from 2 KeV down to a few volts, elastic 
collisions will contribute little to energy loss (see below) . 

2.2.2. Energy Losses Below 20 eV 

Equation (5) will describe the energy versus time down to a 
few eV. The last inelastic collision will take place at the first 
excited state of He which is 19.6 eV. The impacting electron 
will have energy somewhat over 20 eV, and will possess w 1 eV 
after the collision. At 19.6 eV the inelastic cross section 
becomes zero, and below this energy the electrons can only lose 
energy by elastic collisions. The elastic collisions will be of 
two kinds — against neutrals (assumed at around room temperture) 
and against other electrons or icms, where long range Coulomb 
type collisions take place. 

2.2.3. Elastic Collisions with Neutrals 

The rate of change of electron energy versus time for one 
electron from collisions with neutrals is 


6 


d£ 

dt 


-K e V 


(7) 


e 


where K is the fraction of its energy lost in the collision, 
and on the average 


K = 2ra/M = 1/3000 

for electrons in He . The collision frequency here is the 

elastic collision frequency and is equal to Prom equations 

(7) and (3) an energy versus time relation is obtained when 
elastic collisions with neutrals dominate, 



and if the electron starts with an energy ei then 



( 8 ) 


(9) 


Here ^ is a characteristic time to reduce the energy to ei/4: 



( 10 ) 


Assuming Zi ^ 1 eV and X , the mean free path for electron-neutral 
collisions in He^ at 1 torr is 17.6 x 10”^ (ref. 3), then = 

2.3 X 10“^ m at one atmosphere, and 

T 55S 10” ' sec 
e , n 

for the range 1 to 1/4 eV. The function £(t) approaches £ = 0 

asymptotically. The time required to reduce £ from = 1 eV 

down to 1/40 eV is 24 nanoseconds. Note T^ „ v”\ and the 

0 ^ n j. 

range of £i must be taken into account. 
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2.2.4. Elastic Collisions with Cold Electrons 


A measure of the rate of energy loss of electrons of temperature 

T to electrons of temperature Tj can be obtained from the equi- 

partition time t . Strictly speaking this time is a measure 

© " © 

of the rate at which equipartition of energy is established between 
two groups of particles, namely test particles in e Maxwellian 
distribution and at a temperature T, and field particles at a 
temperature Ti . The collisions are assumed to be long range 
Coulomb collisions. The equipartition time is defined by 

the equation (ref. 4) 



and 


eq 


= 3.3 X 1Q5 g3/2 

n In A 


( sec) ; Ti -J" 0 


( 12 ) 


with e = T (eV), electron density, n in cm"^, and A is the 

ratio of the Debye length to the impact parameter. Assume here 

that Ti = 1/40 eV, << T, and let T be 10 eV and 1 eV respectively. 

The values of In A change little with n and T, and the values 

for T are shown in the table below, 
eq 


Ne cm ^ 

In A 

Values 
'Ti = 

of T 

eq 

1 eV 

in sec for 
Ti = 10 eV 

o 

H 

10.7 

3.0 

X 

10" 

10“ ® 

10 

9.6 

3.4 

X 

10"® 

10~® 

10^ ^ 

8.5 

3.8 

X 

10- ® 

10"^ 

lO'"* 

7.3 

4.4 

X 

10"“ 

1.4 X 10“ ® 


S 



For n = 10^^ cm ^ , Ti =1 eV, the time for the hot electrons 

to be cooled by electrons at room temperature is of order 

For energies less than 1 eV the electron-electron collisions 

dominate as t « . 

eq r 

2.2.5. Elastic Collisions with Ions 

The ions, being of equal mass with the neutrals, should be 

quickly cooled, and would then cool the electrons. However, the 

equipartition time for the electrons to be cooled by ions at room 

temperature, t . is far greater than t for electron- 

^ ' eq,e-i eq,e-e 

electron cooling (ref. 4): 


T /x 

eq,e-e' eq,e-x 


ra/M^ = 1/6000 


(13) 


where is the mass of the ion, and the cooling of electrons 

by ions can be neglected. 


2.3. Loss Mechanisms and Plasma Containment Time 

Experimentally, the pulses are a «l/2 ms long, but the 

electron containment time, x could be must shorter. Rough 

c ... 

estimates of x can be made by (1) assuming the electrons 

c> 

diffuse out of the plasma to the walls, or (2) by assuming that 
recombination dominates the life time. 

2.3.1. Containment Time by Diffusion to the Walls 

The containment time, assuming diffusion of electrons to the 
wall, is given by 


where x is a characteristic distance w 1 cm, and D„ the 

a, 

ambipolar diffusion coefficient. The diffusion coefficient for 
helium has been measured (ref. 5) as - 700 ± 50 cm^ sec"*^ at 

1 torr, so X « 1 sec at 760 torr, much longer than the observed 

C • ' 

pulse length. Hence, diffusion is not important. 
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2.3.2. Containment Time by Recombination 

The recombination rate depends on whether He"*" or Het dominates. 
For the former, three recombination processes can occur, namely 

(1) radiative, (2) collisional radiative (He + e + e ->■ He + e), 

+ * 

or (3) neutral stabilized (He •. e + He He + He). Processes 

(2) and (3) are orders of magnitude faster than (1). However, 
with the type of discharges here it seems He! should dominate. 

In addition to processes analogous to (2) and (3), dissociative 
recombination (Hea + e -j- He + He) is also possible. The 
recombination rates of He from (2) and (3), and the recombina- 
tion of Hea , turn out to be the same order of magnitude, which 
yields electrons in ^ terms of n the electron 

density. 


2.3.3. Steady State Density of Electrons 

In the steady staj:e, the production rate of electrons 

■ (*).... 

and 


= 2 X iq2o el/cm~Vsec"^ 

V^Vprod 

from section 1.1. The solution of equation (14) yields 
n « 10 ^ ^ cm” ^ 

and the containment time, t , defined by 
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is of order 10“’ sec, and much less than the pulse duration of 
5 X 10“"* sec. The values of n and t are roughly the same 
irrespective of whether He or He 2 dominates. 

2.3.4. Neutral Gas Temperature 

The final electron temperature, T, should approach the 
temperature of the neutral gas, which is now estimated. Assume 
as above the neutron flux from the moderator is 10^’ per cm^ per 
sec, their average range is 6 cm at one atmosphere (ref. 1), and 
each neutron delivers 0.76 MeV . On the average, 34 eV are 
required per ionization event, and a further 8 eV goes into 
kinetic energy of the electrons. Of the 34 eV, 24 eV goes into 
ionizing, which is recovered after recombination, and probably 
ends up as metastable energy which provides lasing. The remaining 
10 eV goes into excitation and both these contributions end up as 
photons. So, only 8 eV, or 20 percent of the (34 + 8) eV, ends 
in electron kinetic energy, which in turn heats the neutral gas 
(section 2.2.3.). The rate of kinetic energy delivery to the 
neutral gas is therefore 

10^’ X 0,76 X 10® X 0.2/6 = 2.5 x 10^^ eV per cc per sec. 

The pulse length is 5 x 10“"* sec so the kinetic energy per cc is 
1.3 X 10^® eV. In one cc, the number of neutrals is 2.7 x 10^®, 
so each atom gains 0.048 eV, equivalent to a temperature rise 
~ 500 K. 

Again, for an average pulse, the amount of power absorbed by 
one cc of gas per sec « 1 KW (ref. 1). With a pulse duration of 
5 X lO”"* sec, 0.5 joules/cc is deposited, and if 20 percent goes 
into kinetic energy, the temperature of the neutrals should be 
raised by 150 K. The numbers 500 and 150 K are in order of 
magnitude only. 
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2.4. Electron Energy Distribution Function 

2.4.1, Overall Description of the Average Energy of One 
Electron Versus Time 

Consider one electron created at a kinetic energy « 1 KeV. 

It energy first decays rapidly by inelastic collisions with a 
characteristic time and by elastic collisions with a 

characteristic time T „ (T. T „ « In addition, cold 

6 j XX HX € f XX X 

electrons can cool the hot ones. Strong electron-electron cooling 

does not take place until the energies are as low as “ v^. 

However, at £ « 1 eV, t w T , although there are 10® 

eq,e-e e,n 

times more neutrals than electrons. The characteristic times 
for cooling are much less than the plasma containment time. 

The electrons can gain energy from other electrons when at 
low energies by Coulomb collisions (section 2.2.4.), and by 
collisions of the second kind with excited neutrals (superelastic 
collisions) . However, the density of excited neutrals is probably 
of the order of the electron density and much less than that of 
unexcited neutrals; so collisions of the second kind are neglected. 

The picture then is the electrons start off in a distribution 

such as figure 1, and then continuously lose energy until at low 

energies. For energies below 1 eV and n w 10^®, t „ « 

j ©“S c 

and therefore at low energies a Maxwellian distribution of elec- 
trons is attained, with a temperature approximately that of the 
neutral gas. Note that only electron-electron collisions (t o_o) 
lead to a Maxwellian distribution, while electron-neutral elastic 

collisions (T ) need not. The latter reduce the electron energy, 
e , n 

the former redistribute it. This model enables an estimate of 
the steady state electron energy distribution to be made. 

2.4.2. Simple Intuitive Description of the Steady State Electron 
Energy Distribution Function 

The intent of this section is to provide simple intuitive 
arguments to indicate the essential features of the electron 
energy distribution function and compare them with that of Hassan 
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(ref. 9). The duration of a pulse is 5 x 10“**, and the lifetime 
of an electron w 10“^ sec; hence we assume there is a continuous 
constant source of electrons for the duration of the pulse. Each 
electron then starts at a high energy, loses energy with time, 
becomes Maxwellian at low energies, and is recombined. 

The density, dn, of electrons in the energy interval e to 
e + de is proportional to the amount of time, At, each electron 
spends in this interval, 

dn = At 

Now for one electron, 

limit ^ _ / de \ 

Ae^O Ae \dt/ 


where is the rate of change of energy versus time for an 

average electron. Therefore, 


dn = f(e)de 



where f(e) is the steady state electron energy distribution 
function, and 


f (e) « 




We now construct f(e) versus e on a logarithmic plot (fig. 2). 
Starting from high energy, the region DE shows a sharp drop 
for increasing e, the initial distribution of electrons is 
the sum of curves similar to figure 1; the exact shape of DE 
is not discussed further. Below about 35 eV, the electron energy 
loss is mainly from inelastic collisions, where de/dt from 
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equation (4); hence in CD, f(e) The lower limit of 

the region C is around a volt (see section 2.2.2.). Near C 

there is less and less contribution from excitation, and in BC 

electron-neutral collisions dominate. Here equation (8) indicates 

de/dt ® so f(e) « Below B electron-electron 

collisions dominate, and as t << t , a Maxwellian distribution 

eq c 

is formed, with its peak at the temperature of the neutral gas 
at some hundreds of degrees Kelvin. 


The values of the ordinate in figure 2 are not absolute, 
and an order of magnitude calculation enables a comparison to be 
made of the average density, electrons per unit energy 

CD. Assume AB extends over a range 
CD over a range of = 30 eV. Then, 


in AB versus n^^ in 


= 1 eV, and 


^AB _ ^CD 

^ AB ~ t^T^ ^^CD 


'AB 


'CD 


where the times t 

spends in AB and CD. 

tion time, and t^-p, = T.^ 

CD in 


CD 
Now 


correspond to those an average electron 


AB "'c = 10 
= 0.4 X 10"'“ sec. 


- 7 


sec, the recombina- 
Theref ore 


n 

n 


CD 


10^ or lO"* 


The relative scale on the ordinate of figure 2 conforms with this 
ratio . 

The essential details fo figure 2 agree with an electron energy 
distribution curve of Hassan (ref. 6) derived from the Boltzmann 
equation (fig. 3). Hassan actually plotted f(e)/e^/^, and the 
curve for f(e) is due to Meador and Weaver (ref. 7). The points 
B and D in figure 2 occur in the same energy range as figure 3. 
Between B and D Hassan 's curve shows f(e) ^ whereas 

figure 2 shows f(e) « near D, and « near B. Both 

curves show a Max^vellian distribution to the left of B with 
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about 10,000 times more electrons in the Maxwellian distribution 
below 1 eV than there are in the non-Maxwellian distribution 
above it . 


3. FUTURE WORK 

It would be desirable to measure the electron density and 
temperature, and a study has been made on the feasibility of using 
probes. A paper on high pressure probes by Cozens and von Engel 
(ref. 8) implies that the measurements are possible. Some pre- 
liminary experimental tests are to be performed in collaboration 
with N. Jalufka on an electrically created plasma, and if the 
probes work well, possibly the values checked using other nonprobe 
methods . 

The work shown in section 2 ±'j preliminary and refers to the 
electrons only. Future studies will include mathematical models 
which describe the ion processes including laser action. 
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